A band structure scenario for the giant spin-orbit splitting observed at the Bi/Si(lll) interface 
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The Bi/Si(lll) (\/3 x \/3)R30° trimer phase offers a prime example of a giant spin-orbit splitting of the 
electronic states at the interface with a semiconducting substrate. We have performed a detailed angle-resolved 
photoemission (ARPES) study to clarify the complex topology of the hybrid interface bands. The analysis of 
the ARPES data, guided by a model tight-binding calculation, reveals a previously unexplored mechanism at 
the origin of the giant spin-orbit splitting, which relies primarily on the underlying band structure. We anticipate 
that other similar interfaces characterized by trimer structures could also exhibit a large effect. 



I. INTRODUCTION 

The normal spin degeneracy of the electronic states of non- 
magnetic solids is lifted by the spin-orbit (SO) interaction in 
crystals lacking an inversion center (Dresselhaus effect) |IT]|2l. 
A similar effect was predicted theoretically by Rashba and 
Bychkov (RB) for a two-dimensional electron gas (2DEG) at 
a surface or an interface which exhibits a structural surface 
asymmetry (SSA) |i3J. Although the model was originally 
motivated by semiconductor heterojunctions, split bands were 
first observed by angle-resolved photoemission spectroscopy 
(ARPES) on metal surfaces |4-10|. Like its atomic counter- 
part, the RB effect has a relativistic origin, namely the cou- 
pling of the spin to the magnetic field which appears in the 
rest frame of the electron. In the free-electron limit consid- 
ered by RB, the parabolic dispersion is split, as in Fig. 1, into 
two branches of opposite spin: 

E^{k)^—±aRk, (1) 

where k is the magnitude of the electron momentum in the 
plane of the surface, and m* the effective mass. The Rashba 
parameter q;r is proportional to the gradient of the surface 
electric potential, and defines the strength of the RB effect. 
The SO-splitting of the two branches can be quantified by 
their momentum offset 2fco = 2aiim* /h^ or, equivalently, 
by the Rashba energy Eb, ~ /i^fc§/(2m*), the difference be- 
tween the band minimum (to* > 0) or maximum (m* < 0) 
and the crossing point of the two branches at A; = 0. In a 
more realistic approach, the band splitting depends not only 
on the surface potential gradient, but also on the atomic SO 
parameter, and on the asymmetry of the electron wavefunc- 
tions [10-13|. 

Interest in the RB effect has been revamped by observations 
of a giant SO-splitting in surface alloys formed by a high-Z 
element - Bi or Pb - at the Ag(lll) surface OKBl. The 
unusual strength of the effect has prompted a reassessment 
of the various factors contributing to the effect. It has been 
suggested that additional components of the surface potential 
gradient within the surface, reflecting the anisotropic charge 
distribution, are probably important |16|. Independent studies 




FIG. 1 : (color online) Schematics of the RB SO split bands for a 2D 
free electron gas 



have stressed structural aspects, namely the role of relaxation 
and buckling of the topmost layer in defining the hybrid states 
El. 

Metallic surface alloys with a giant SO splitting are po- 
tentially interesting for spintronics applications. The present 
challenge is to make them compatible with semiconductor 
technology |T8l,T9 |. There has been encouraging progress in 
this direction, and several attempts have been made to grow, 
on Si(lll) substrates, thin layers which support SO-split 
bands at their surface |'20'-'24l. Recently, a giant spin-splitting 
with no buffer layers was demonstrated for the isostructural 
Bi/Si(lll) 12511^ and Bi/Ge(IlI) W\ interfaces. In the 
same line, metallic spin-split surface states were observed in 
the related Pb/Ge(ll 1) system |28 1. In all cases, the electronic 
structure of these interfaces is more complex than that pre- 
dicted by the simple RB model. Although first-principles cal- 
culations reproduce the experimental results, they suffer from 
a certain lack of transparency. This leaves room for a simpler 
but more direct approach which can help in the interpretation 
of the experimental data, and thus contribute to clarify the un- 
conventional properties of the electronic states. This was the 
motivation of the present work, which compares the results 
of a detailed experimental band mapping of the Bi/Si(lll) in- 
terface by ARPES, with simple models of the band structure 
in the presence of a RB-like interaction. In particular, a para- 
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FIG. 2: (color online) The structure of (a) the monomer (a—) and (b) 
the trimer (/? - i/3 x \/3)R30° Bi/Si(l 1 1) phases. The size of the Si 
atoms indicates their distance from the surface. Solid (dashed) lines 
indicate the 1x1 (\/3 x \/3) primitive unit cell. The horizontal and 
vertical black lines follow the mirror planes of the overlayers. (c) The 
1x1 (large hexagon) and (\/3 x \/3)R30° (small hexagons) SBZs. 
The circles are constant-energy contours for RB paraboloids centered 
at each equivalent F point of the latter, and the arrows indicate the 
spin polarization. The reciprocal lattice vectors Gi, G2 and G3 = 
Gi — G2 are those considered in the NFE model of subsections III 
(B) and (C) 

metric tight-binding scheme provides a satisfactory qualitative 
description of the data, and suggests a possible new mech- 
anism to achieve a large spin polarization, which is closely 
connected with a characteristic feature of the band structure 
of the interface. 



II. EXPERIMENTAL DETAILS 

The Si(l 11) substrate (Sb-doped, resistivity O.Olfi.cm) was 
flashed at 1200°C by direct current injection, and then cooled 
slowly in order to obtain a sharp low-energy electron diffrac- 
tion (LEED) (7 X 7) pattern. The (^3 x vl)R30° Bi/Si(l 11) 
interface was prepared by deposition of 1 monolayer (ML) of 
Bi from a calibrated EFM3 Omicron source on the the sub- 
strate at RT followed by a mild annealing. ARPES spectra 
were acquired at 70 K and 21.2 eV photon energy, with a 
PHOIBOS 150 SPECS Analyzer equipped with a monochro- 
matized GammaData VUV 5000 high brightness source. The 
Fermi level position was determined from the Fermi edge of a 
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FIG. 3: (color online): 2"'' derivative ARPES intensity map along 
the rMr direction for the monomer (a) and the trimer (b) phases of 
Bi/Si(lll). 

polycrystalline Au sample. 

III. RESULTS AND MODEL CALCULATIONS 

A. ARPES Measurements 

The Bi/Si(lll) interface exhibits two different structures 
with the same (^3 x V3)R30° Bi/Si(lll) periodicity: a 
monomer structure (a— phase) for a coverage of 1/3 ML, and 
a trimer structure (/3— phase) for 1 ML coverage. They are 
illustrated in Fig. 2 (a) and (b). According to the widely ac- 
cepted T4 model 1 29 1, both monomers and trimers are cen- 
tered above the 2"*^ layer Si substrate atoms. Figure 2(c) 
shows the surface Brillouin zones (SBZ) of the unrecon- 
structed Si(lll) surface, and of the (V^ x V3)R30° super- 
structure. In the following we will always refer to the latter 
The a and /3 phases have quite different band structures, and 
can be easily distinguished. This is illustrated in Fig. 3, where 
(second derivative) ARPES intensity maps of the two phases 
are compared around the M point. In agreement with previ- 
ous studies, the a— phase shows a rather flat surface state with 
a broad minimum at M |30|, while the /?— phase shows two 
symmetrically dispersing features crossing at M |3T!l32l. In 
the rest of the paper we only consider the trimer /?— phase, 
which is the most interesting in the present context. The sys- 
tem has a threefold rotation axis, and three mirror planes per- 
pendicular to the surface. One of these mirror planes is par- 
allel to the [112] direction, and also to the FKM direction of 
the SBZ {ky axis in Fig. 2 (c)). The other mirror planes are 
rotated by 120° around the z axis. The overlayer symmetry 
is identical to the one of the substrate and corresponds to the 
plane group p31to. The characteristic band crossing at M is 
the signature of a peculiar RB-type SO-splitting with a mo- 
mentum offset fco of 0.126 A^^ and a Rashba energy of 140 
meV 1251 126< . The large splitting has been previously asso- 
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FIG. 4: (color online) (a) to (c) Second-derivative ARPES intensity 
maps showing the band structure along the fc-space cuts (a) to (c) 
of panel (d). The dashed curves and arrows are guides to the eye, 
inspired by results obtained in Refs. 125 1 1261 . and highlight the dis- 
persion of the Si split branches. The thick solid line marks the edge 
of the projected Si bulk gap. Images (a) and (b) intersect image (c) 
respectively along the vertical dashed and solid lines. 



ciated with the inversion asymmetry induced by the trimers. 
We will see later that a somewhat different interpretation is 
possible. 

The ARPES intensity maps of Fig. 4 illustrate the disper- 
sion of three surface states - labelled Si to S3 as in Ref. |26] 
-. All three states are predicted to be spin-polarized (26\. Si 
exhibits a large splitting around M, as already shown by Fig. 
3, and a peculiar anisotropic dispersion around that point. The 
experimental dispersion along the TMT line (Fig. 4 (a)) shows 
a hint of the two branches 1251 |26l . which split away from the 
crossing point at M. The dispersion is highlighted by dashed 
guides to the eye, which are consistent with the results of Syn- 
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FIG. 5: (color online) (a) Experimental CE maps around M. The 
corresponding binding energies are 0.53, 0.68, 0.84, and 1.02 eV 
for panels (I) to (4). They are indicated by horizontal dashed lines 
in the intensity maps of panels (b) and (c), taken along the KMK and 
FMr directions. In the gray-scale plots, highest intensity is black. 

chrotron Radiation (SR) studies, where the individual Si com- 
ponents could be more clearly resolved along FMr ll26l . The 
two S 1 branches merge again approaching the T point, where 
they cannot be resolved from S2. The energy splitting of the 
two branches increases away from the FMF high-symmetry 
line, as shown by Fig. 4 (b) which shows the dispersion along 
the parallel cut (b) of Fig. 4 (d) (i.e. ky 7^ 0). In support of 
our previous discussion on the experimental results of Figs. 4 
(a) and (b), the Si split branches are easily identified in Fig. 
4 (c), which shows a cut parallel to the KMK line. Along this 
line the intensity of S2 is very small due to ARPES matrix 
elements. The images of Fig. 4 (a) and (b) intersect the per- 
pendicular cut (c) along respectively the vertical black dashed 
and the red solid lines. 

Figure 5 illustrates the unusual topology of surface band Si 
around the M point, in a region of fc-space where it is well 
separated from other surface or bulk-derived features. Panels 
1 to 4 of Fig. 5 (a) show constant-energy (CE) cuts taken at 
increasing binding energies between 0.53 and 1.02 eV, corre- 
sponding to the horizontal dashed lines in the intensity maps 
along the KMK and FMF directions of panels (b) and (c). 
Starting at the highest binding energy, the CE maps show two 
intersections along KMK, symmetrically located with respect 
to M. They get closer at lower binding energy, following the 
negative-mass dispersion of Fig. 5 (b). The two intersections 
finally merge at M for Eb ^ 0.68 eV. This is the maximum 
of the dispersion along KMK, and the crossing point of the 
two SO-split branches along the perpendicular FMF direc- 
tion. Panel 1, taken above this energy, indeed shows non- 
intersecting CE contours. 



4 



The data of Fig. 5 reveals that the topology of Si is quite 
different from that of Fig. 1, predicted by the usual Rashba 
model for a free-electron band centered at the T point. The 
degeneracy of the SO-split branches at M is required by a 
combination of time-reversal and translational symmetry. On 
the other hand, the line of (near) spin degeneracy of Fig. 5 (b) 
finds no correspondence in the simple Rashba model. The ex- 
perimentally observed large difference of the energy splitting 
along the two high-symmetry directions is well captured by 
first-principles calculations Il25ll26i . 

B. An Isotropic Nearly-Free Electron Model 

We will now attempt a comparison of the experimental data 
using simple theoretical models, to gain further insight in the 
unusual dispersion of the SO-split bands. The minimal re- 
quirement for any model is that it should include both the 
Rashba-like interaction and translational invariance. The sim- 
plest approach satisfying this condition is an isotropic nearly- 
free electron (NFE) model. This model is schematically il- 
lustrated in Fig. 2 (c). RB paraboloids (m* > 0) are cen- 
tered at all equivalent F points, and the CE lines are concen- 
tric circles representing the inner and outer SO-split branches. 
The spins exhibit a vortical structure around the F points. 
Two paraboloids centered at two adjacent SBZ intersect in the 
common KMK Brillouin zone boundary along two parabolas, 
one at lower energy for the outer SO branch, and a second at 
higher energy for the inner 



The Rashba hamiltonian for a free-electron is (3): 

H^B{k) ^ an{(T y. k)^ , (2) 

where cr is the vector of Pauli matrices. In a representation 
where the basis states are |fe f), \k \), and spin projections 
refer to the z axis, the corresponding matrix is: 

(h'^k^/2m aii{ky + ik^) \ 
■ (3) 
ctB,{ky — ikx) h^k"^ /2m J 

Diagonalization of this hamiltonian generates SO-split 
paraboloids centered at F as in Fig. 1. The corresponding 
'spin-up' and 'spin-down' eigenstates refer to a quantization 
axis ee = X k//c, which is always perpendicular to k, 
i.e. tangential to the constant-energy circles. These states are 
therefore 100% in-plane polarized, with a purely tangential 
spin polarization P - opposite on the two branches - rotat- 
ing around F. This is easily generalized to include the lat- 
tice periodicity. Since we are mainly interested in the band 
structure near E p around the M point, it is a good approxima- 
tion, to consider only the first SBZ and three adjacent zones 
centered at the lattice vectors Gi = (1,0), G2 = (0,1) 
and G3 = Gi — G2, as in Fig. 2 (c). The basis vectors, 
again referred to the z axis, are \k X), \k \), |(fe + Gi) '\), 
\{k + Gi) i), |(fc + G2) t), \{k + G2) ;), \{k + G3) t), 
I (fc + G3) I) The truncated NFE hamiltonian matrix is then: 



/i?RB(fe) K)l 1^02 1^03 \ 

^0*1 i^RB(fc + Gi) 1/12 1^13 

1^0*2 ^1*2 irRB(fc + G2) V23 

\ ^03 Viz ^^23 HRB{k + Gs)J 



The diagonal 2x2 building blocks now generate SO-split 
paraboloids centered at F, and at Gi, G2 and G3. The off- 
diagonal 2x2 blocks describe the interaction between states of 
equal spin on the various paraboloids ll33l . The hybridization 
strength is, as usual, the corresponding Fourier component of 
the crystal potential V{r) defined by: 

V{r)=Y.^G,e'''"^ . (5) 

i 

For instance: 

It is easy to show that (6) is equivalent to an interaction of the 
form V^(k, k + Gi) = VgiCos((5/2), where S is the angle be- 



tween the two polarization vectors. In our case all hybridiza- 
tion terms are equal: Voi — Vg, and Hj^fe contains only 
the two parameters and Vq. The binding energy of the 
paraboloids at F, or equivalently the Fermi level position, are 
adjusted to fit the experimental data of Fig. 4. 

Figure 6 (a) illustrates the predictions of the NFE model 
along the FMF direction, in the limit \Vg\ — 0. The outer 
branches cross at the M point, with opposite spin polarization. 
Along the perpendicular KMK direction (Fig. 6 (c)), their 
intersection is a parabola dispersing upwards from M. The 
inner branches of the paraboloids similarly cross above Ep. 
The model yields a constant momentum separation between 
the SO-split branches of each paraboloid, and does not capture 
the experimentally observed fc-dependent splitting. Moreover, 
the sign of the dispersion along the SBZ boundary is opposite 
to that of Fig. 5 (b), and the opposite spin states are strictly 
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FIG. 6: (color online) (a), (b) Band dispersion of the isotropic NFE 
model along the rMT for |Vg| = and |Vg| = 0.3 eV. Arrows 
indicate the opposite (in-plane) spin polarization of the two branches, 
(c), (d) Same for the KMK high-symmetry direction. The two spin 
states are degenerate along the SBZ boundary for |Vg| =0. Notice 
the different scales in (a), (b) and (c), (d). 



degenerate, for any value of 

Figure 6 (b) shows that the main effect of a finite lattice po- 
tential is the opening of energy gaps at the crossing of bands 
with parallel spins, i.e. at the crossing of the outer branch of 
one paraboloid with the inner branch of the paraboloid cen- 
tered at an adjacent SBZ. No gap opens when either the outer 
or the inner branches cross at the M point, because their spins 
are opposite there. This is consistent with the requirements of 
time-reversal symmetry. At the other crossing points along the 
SBZ boundary the spins are not strictly opposite, but they are 
nonetheless rather antiparallel, and the hybridization is there- 
fore small. Therefore, as shown in Fig. 6 (d), the two spin 
states are not degenerate as in the case IVgI =0, but their en- 
ergy separation is small along this direction. The model again 
yields a positive effective mass along KMK, in contrast with 
the experiment. The momentum splitting along KMK scales 
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FIG. 7: (color online) (a), (b) Band dispersion of the anisotropic 
NFE model along the KMK SBZ boundary for | Vg | = and \Vg\ = 
0.3 eV. (c) Constant energy contours for a binding energy of 1.1 eV 
showing the anisotropic shape of the spin-split states. Green (light 
gray) and brown (dark gray) colors indicate positive and negative 
values of tangential spin polarization. 



with IVgI and a^l"^, confirming the unconventional character 
of the underlying mechanism. It should also be noted that al- 
though the M spin degeneracy is fundamental, the degeneracy 
predicted at K is accidental. It is lifted when further reciprocal 
lattice points are included. 



C. An Anisotropic Nearly-Free Electron Model 

The positive effective mass along the Bragg planes pre- 
dicted in Fig. 6 (c), (d) is a direct consequence of the sim- 
ple circular CE contours of the isotropic NFE model. More 
complex contours are however possible in anisotropic 2D sys- 
tems. Blossom-like contours have been predicted for Rashba 
systems |14, 16 1, and experimentally observed for the BiAg2 
||25| and SbAg2 I34l surface alloys. Concave CE contours 
are a generic effect of an in-plane anisotropy of the potential 
and they are not limited to surface alloys. For example, a 2D 
Dirac fermion state at the surface of the topological insula- 
tor Bi2Te3 has been shown to exhibit a snow-flake-like Fermi 
surface [35 1. 

k • p theory has been used to calculate higher-order terms 
in the effective Hamiltonian of a topological insulator with 
i?3m symmetry |35 1. Following these results we introduce an 
anisotropy in H^PEik) as: 



(7) 
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FIG. 8: (color online) Constant energy maps as measured by ARPES (middle panel) and according to the predictions of the isotropic (left 
panels) and the anisotropic (right panels) NFE models described in the text. Both models are presented for | Vg] — and \ Vg\ = 0.3 eV. 
Green (light gray) and brown (dark gray) colors indicate positive and negative values of tangential spin polarization. The arrows are a sketch 
of the predicted in-plane projection of the spin polarization around the M point. 



with 

H'{k) = ^{{ky + ik^f + [ky - ik^f)a, . (8) 
This yields: 

E^{k)^ ± \/("i?)^fc^ + c2fc6cos2(36') , (9) 

2m* 

where c is an anisotropy parameter, and is the in-plane angle 
from the FM direction. For small values of c this expression 
reduces to the free-electron case with RB splitting (Eq. (1)). 
The in-plane asymmetry could be self-consistently included 
starting from Eq. (5), but Eq. (7) provides a minimal alter- 
native model with the single parameter c. The resulting band 
dispersion agrees well that of the anisotropic 2DEG proposed 
by Premper et al. II16J . It is shown in Fig. 7 (a) for \ Vg\ — 0, 
and (b) for = 0.3 eV. The coiTesponding parameters are 
summarized in the Appendix. The model correctly predicts a 
negative effective mass along KMK. Again, the crystal poten- 
tial induces a weak RB splitting. The dispersion along FMF 
is essentially identical to the isotropic case and, as the latter, 
it overestimates the spin-splitting of the bands, especially far 
from M (not shown). The CE contours evolve continuously 
with increasing energy from circular to blossom-like, as in 
Fig. 7 (c). The sixfold symmetry is the result of the threefold 
rotational symmetry and time-reversal symmetry. It should 
be noted that the spin polarization has only an in-plane tan- 
gential component for an isotropic 2DEG, whereas a sizeable 
out-of-plane and a small radial component are present for the 
anisotropic case |16|. Fig. 8 compares the experimental CE 
contours around M with the predictions of the isotropic (left 



panels) and the anisotropic (right panels) NFE models. The 
latter describes reasonably well the data for \ Vq\ —0 but fails 
to reproduce the dispersion (Fig. 6 (a)) due to the overesti- 
mation of the momentum splitting at fc-points far from M. A 
finite lattice potential does open an energy gap, but it perturbs 
the band structure and the agreement is completely spoiled. 

A closer examination of the CE contours reveals another 
subtle inaccuracy of the anisotropic NFE model. The spin po- 
larization symmetry, determined by the mirror plane of the 
trimer configuration 1. 361 , is not in agreement with the tip ori- 
entation of the outer CE contour of Fig. 7 (c). This inconsis- 
tency is removed by the tight-binding model considered in the 
following subsection. 

D. An Empirical Tight-Binding Model 

The covalent character of the bonds and the semiconducting 
nature of the Bi-Si(lll) interface suggest that local orbitals 
may be a better starting point. We have performed an empir- 
ical tight binding (TB) calculation which is able to reproduce 
the main experimental features. In order to limit the complex- 
ity of the calculation, the model considers a single orbital per 
atomic site with sp^ symmetry. While this is an approxima- 
tion, we expect contributions from other orbital symmetries to 
be small in the energy range of interest. This is supported by 
recent results for the isostmctural Bi-Ge(lll) interface |27|. 
The primitive unit cell contains three Bi atoms, labelled a, b 
and c in Fig. 9. In the same figure the five inequivalent hop- 
ping terms are indicated by arrows. All the other terms can 
be generated by symmetry. The Si(lll) substrate is only in- 
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FIG. 9: (color online) Schematics of the cell used for the TB calcu- 
lation, with the definition of the inequivalent transfer integrals. 



directly taken into account through the effective hopping pa- 
rameters. A calculation of the transfer integrals is a non trivial 
computational task, which clearly goes beyond the scope of 
this work. Therefore we defined them in a purely phenomeno- 
logical way, assuming an inverse power-law dependence of 
the distance d between two centers: V{d) = ad~^. The pref- 
actor a determines the bandwidths, while the exponent b deter- 
mines details of the dispersion. There is obviously no angular 
dependence for spz states. The Bi-Bi distance within a single 
trimer was set to 2.6 A, which is very close to literature val- 
ues Il25ll37]| . In the actual calculation we included interactions 
up to 4* nearest neighbors. The required overall resemblance 
with the experimental dispersion significantly limits the ac- 
ceptable parameter space. The chosen values are summarized 
in Table II of the Appendix. 

Figure 10 illustrates the results of the TB model before the 
inclusion of a RB interaction. The calculation yields three 
bands, corresponding to the three orbitals per unit cell. The 
two higher lying states are shown for the FMP (a) and FKM 
(b) directions in Fig. 10. They can be associated to the ex- 
perimental features of Fig. 4. The model correctly predicts 
a double degeneracy at the F and K points, independent of 
the parameter values. A double degeneracy is imposed by the 
Csv symmetry of these points for the trimer structure. The M 
point has a lower symmetry (Cih) and therefore no degener- 
acy is expected there. An interesting result is the presence of 
two maxima for Si on both sides of M along FMF. This hall- 
mark of the trimer structure, which is not observed in simple 
hexagonal structures, has been repotted for other similar sys- 
tems [38', "39"!. We shall see below that it plays a potentially 
important role in the appearance of a giant SO splitting in 
Bi-Si(lll). 

We consider next the effect of the SO interaction, by adding 
a RB term to the TB hamiltonian Ii40h43i : 

HjB = ^y4cj^+« ^y4('^xdij)^Cj,' (10) 

<i.,j> <i,j>s,s' 
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FIG. 10: Calculated band dispersion of the two higher lying states 
according to the TB model without RB interaction, along the FMF 
(a) and the FKM (b) high-symmetry directions. 

The first term is the usual spin-independent TB hamiltonian, 
while the second term is the appropriate TB form for the 
Rashba interaction. cl^{cis) is the creation (annihilation) op- 
erator of an electron with spin s (^ or I) on atomic site i, 
Vij and Xij are the transfer integrals and the SO coefficients. 
The latter are generated by a similar power-law function of the 
distance, but with independent parameters (see Table II of the 
Appendix), cr is the vector of the Pauli matrices and dy is the 
vector connecting site j to site i. 

The six basis vectors of a site-spin representation are, with 
obvious notation, \a t), \b t), |c t), \a |), \b |), |c |). In this 
representation, the hamiltonian matrix has the form: 

Hf) is the 3x3 spin-independent TB hamiltonian which de- 
scribes states of equal spin. Hr is a 3 x 3 matrix generated 
by the second term of Hjs, which describes the coupling of 
electrons with opposite spins. Our method is essentially akin 
to the TB model of Ref. [TT]and can qualitatively describe the 
SO split bands, namely their topology. 

Figure 1 1 shows the band structure of the TB model with 
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FIG. 11: (color online) Calculated bands for the TB model with SO 
(see Table II of the Appendix for the parameters) for: (a) the FMF 
dire ction; (b) the parallel cut b of Fig. 4 (d) (ky = 0.05 A"^) (b); (c) 
the FKM direction. Arrows denote the main component of the spin 
polarization of the SO-split branches. The polarization difference is 
not 100% due to additional radial and out-of-plane components. 

RB interaction, for the parameters which best reproduce the 
experimental results [44 1. The first obvious effect of the RB 
interaction is that all states are now split. Both Si and S2 
are doubly-degenerate at F, and exhibit an isotropic disper- 
sion around this point (Fig. 1 1 (a)). By moving out of the 
high-symmetry line (Fig. 1 1 (b)) the degeneracy is lifted, as 
expected for the usual RB scenario of Fig. 1. The spin-split 
branches cross again at M, as required by time-reversal sym- 
metry. Here, the dispersion of Si is strongly anisotropic. The 
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FIG. 12: (color online) (a) CE ARPES contours. The dashed curves 
on the experimental maps are guides to the eye. Green (light gray) 
and brown (dark gray) colors indicate positive and negative values of 
tangential spin polarization, (b) CE contours for the TB model. The 
arrows indicate the in-plane projection of the spin polarization. The 
energies were adjusted to correspond to the experimental values. 



splitting of Si along the FMF direction is much larger than 
that of S2. It is also much larger than the splitting of Si around 
F. By contrast, it is small along the SBZ boundary KMK (Fig. 
1 1 (c)). All these features of the band structure agree with the 
ARPES results, and also with the results of first-principles cal- 
culations II25I l26l . The calculated energy difference E/j(M) 
between the band crossing at M and the band maximum is 
nonetheless smaller than the experimental one. 

A comparison of Fig. 10 (a) and Fig. 11 (a) shows that 
the large SO splitting of Si is a consequence of the split max- 
ima on opposite sides of M, rather than the result of a large 
SO coupling. In a way, the conditions for a large momen- 
tum splitting are already present in the band structure, and the 
main effect of the SO interaction is to split in energy the two 
subbands. Indeed for S2, which has a maximum at M in the 
parent structure, the Rashba splitting is small. The origin of 
the large momentum separation is therefore rather different 
from that of the 'giant' splitting observed at metallic inter- 
faces lilce BiAg2/Ag(lll) 114|. In the TB frameworlc, it is 
due to the second term of the Hamiltonian (Eq. (10)), i.e. the 
way the trimer arrangement determines the hybridization of 
Bi orbitals with unlike spins. 

Figure 12 illustrates the energy evolution of the calculated 
CE contours near M, across the spin-degeneracy point. The 
experimental contours are reproduced here for a qualitative 
comparison. A fully quantitative comparison is not possible 
due to the already mentioned difference in E;j(M). The en- 
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ergies of the calculated contours were therefore adjusted to 
correspond to the experimental energies of the ARPES con- 
tours. The model yields open contours around M, which are 
in reasonable agreement with the topology of the experimental 
bands. The CE contours are shown on a broader momentum 
range in Fig. 13 (a) and (b) for two energies above, and in 
Fig. 13 (c) for an energy well below the crossing point. The 
shape of the contours is nearly circular near the bottom of the 
band at F (Fig. 13 (c)), and it evolves to a hexagonal and fi- 
nally blossom-like shape at larger energy. This is seen more 
clearly in Fig. 13 (d), where the SO parameters were artifi- 
cially increased to enhance the splitting of the two subbands. 
Interestingly, the blossom-like shape of Fig. 13 (d) is identical 
to the one predicted by the anisotropic NFE model (see sub- 
section III (C) and Ref. 1 16 1 but the tips now correctly point 
along the FMF direction. 



IV. CONCLUSIONS AND OUTLOOK 

We performed a detailed ARPES study of the SO-spUt elec- 
tronic states in the 1 ML trimer phase of the Bi-Si(lll) in- 
terface. We paid special attention to the region of fc-space 
close to the M point, where the topmost Si hybrid surface 
state is both not degenerate with bulk states and distinct from 
other surface-related features. This region is of particular in- 
terest because there Si exhibits a large and non-conventional 
Rashba-like splitting. Energy-dependent constant energy con- 
tours clarify the complex topology of the SO-split states, and 
underline the differences with a standard RB scenario. The 
ARPES data show that the interface has an insulating charac- 
ter, but the Fermi level could be moved into the SO-split bands 
by applying an external electric potential in a back-gated 
structure. It should then be possible to advantageously exploit 
the large momentum separation of the two spin-polarized sub- 
bands in a spin field-effect transistor |fT8l[T9ll . 

We have used the predictions of three simple models for a 
2DEG in the presence of SO interaction as guidelines for the 
interpretation of the experimental results. The comparison of 
the NFE and local-orbital schemes, which proceed from op- 
posite starting points, has a certain didactic value. Moreover 
we were able to assess the limits of the various approaches 
applied to the Bi-Si(lll) case. The ARPES results and CE 
contours are well described by a NFE RB effect in a suffi- 
ciently small region around the F point. Further away from 
r the isotropic NFE model must be refined to include an in- 
plane asymmetry. Near the SZB, the specific symmetry prop- 
erties of the interface determine the characteristics of the SO 
splitting |45|, and only the empirical TB model captures the 
salient features of the electronic structure. The same model 
predicts a peculiar spin texture (Fig. 13 (a)). Hole pockets 
with a non-vortical spin arrangement are reminiscent of the 
teardrop Fermi surface contours of the topological insulator 
Bii xSbx P6l . Within a few meV the pockets develop into two 
connected concentric contours with the same spin polarization 
(Fig. 13 (b)). This new prediction calls for an experimental 




-0.3-o.a-o.i a 0.1 0.2 0.3 -0.3-0.2-0.1 0.1 0.2 0.3 

FIG. 13: (color online) CE contours of the TB model for energies 
above (a), (b) and well below (c) the spin degenerate point at M. 
Arrows indicate the in-plane projection of the spin polarization, (d) 
is the same as (c) for an eightfold increase in the SO parameters. 
Green (light gray) and brown (dark gray) colors indicate positive and 
negative values of tangential spin polarization. 



verification by spin-resolved ARPES. 

Our TB approximation obviously cannot reproduce all the 
details of the complicated electronic structure. More elabo- 
rate TB schemes could be implemented by extending the set 
of local orbitals, by treating in an explicit way the hybridiza- 
tion with the substrate, and by deriving the relevant transfer 
integrals from a direct calculation. However, the actual merit 
of such schemes would be dubious, since the computational 
complexity would approach that of first principles calcula- 
tions, and since the immediate simplicity of the model would 
be lost. 

The main insight from our analysis of the new experi- 
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mental data is the realization that the giant SO splitting at 
this interface is not primarily controlled by a large atomic 
SO parameter, as in the case of BiAg2/Ag(lll) and other 
metallic surface alloys |14|. On the contrary, the effect is 
largely due to a peculiar feature in the band structure, namely 
the presence of symmetrically split maxima around the M 
point. This is a new, unexpected mechanism to achieve large 
spin separation at an interface. The underlying band feature 
is characteristic of the trimer structure, and it has been iden- 
tified at other similar interfaces ll38l [39l . We may therefore 
anticipate that similar large "Rashba"-hke effects could be 
discovered in other systems characterized by moderate SO 
parameters. 

It is a pleasure to acknowledge fruitful discussions with A. 
Baldereschi, J. Henk, I. Rousochatzakis and K. Sakamoto. We 
also thank L. Moreschini and L. Casanellas for help during the 
early stage of this work. E.F. acknowledges the Alexander S. 
Onassis Public Benefit Foundation for the award of a scholar- 
ship. This research was supported by the Swiss NSF and the 
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APPENDIX 

The phenomenological parameters of the models are sum- 
marized in the following tables; 

TABLE I: Parameters of the NFE modes. The potential value refers 
to the 1st Fourier coefficient of the crystal potential V(r). 



parameter isotropic NFE anisotropic NFE 


m* (me) 


0.8 


0.8 


c (eV.A^) 





4.8 


an (eV.A) 


1.1 


1.1 


\Vg\ (eV) 


0.3 


0.3 



TABLE II: Parameters of the TB model. The hopping and SO pa- 
rameters are generated by power law functions of the distance d (i.e. 

ad~''). 









prefactor a 


-2.94 


0.15 (1.20 in Fig. 13 (d)) 


exponent b 


1.13 


0.80 
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